Late Inspiral and Merger of Binary Black Holes in Scalar-Tensor Theories of Gravity 
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Gravitational wave observations will probe non-linear gravitational interactions and thus enable 
strong tests of Einstein's theory of general relativity. We present a numerical relativity study of the 
late inspiral and merger of binary black holes in scalar-tensor theories of gravity. We consider black 
hole binaries in an inhomogeneous scalar field, specifically binaries inside a scalar field bubble, in 
some cases with a potential. We calculate the emission of dipole radiation. We also show how these 
configurations trigger detectable differences between gravitational waves in scalar-tensor gravity and 
the corresponding waves in general relativity. We conclude that, barring an external mechanism to 
induce dynamics in the scalar field, scalar-tensor gravity binary black holes alone are not capable 
of awaking a dormant scalar field, and are thus observationally indistinguishable from their general 
relativistic counterparts. 



PACS numbers: 

Experimental tests in connection with gravitational 
redshift, light deflection, Shapiro time delay and peri- 
helion advance, among others, have increased our confi- 
dence that general relativity (GR) is the correct theory 
of gravity Jj. More compelling evidence is provided by 
binary pulsar observations [2], where the hardening of 
the binary is accounted for in exquisite detail by one of 
the fundamental predictions of GR: gravitational wave 
(GW) emission. Modified gravity theories are nonethe- 
less a possibility [I]. To verify that indeed the podium 
only belongs to Einstein's theory, strong-gravity tests are 
needed, involving for instance core-collapse supernovae or 
the last few orbits and coalescence of compact objects. 
The new astronomy of GW observations will soon deliver 
such opportunities, probing gravity at its strongest grip, 
and thus testing whether Einstein was correct. 

In anticipation of GW observations by LIGO, Virgo 
and other interferometric detectors, exploring what to ex- 
pect from alternative theories is crucial in assisting data 
analysis efforts. Of particular importance is the investi- 
gation of predictions from alternative theories of gravity 
during the generation of GWs in the non-linear regime. 
This calls for theoretical studies only accessible with the 
tools of numerical relativity, which is the focus of the 
present paper. Among the competing alternatives to GR, 
scalar-tensor (ST) gravity [21 E] is one of the most popu- 
lar due to its simplicity, and because of the motivations 
provided by string theory scenarios and explanations to 
dark energy [5] . ST gravity in its simplest form was pro- 
posed about a half century ago and is commonly known 
as Brans-Dicke theory [B]. 

Studies on observational consequences of ST theories 
have mostly focused on compact object binaries in the 
post-Newtonian or extreme-mass-ratio regimes [7l416j. a 
critical prediction of which is the emission of dipole ra- 
diation. The latter depends on the sensitivity, s, of the 



compact objects [T], which is a measure of how suscepti- 
ble the mass of an object is to variations in the local value 
of the gravitational constant (s — 1/2 for black holes 
(BHs), and s < 1/5 for neutron stars |7J. With modi- 
fications to GR entering as the difference As = Si — S2 
of the binary components, binary black holes (BBHs) are 
unaffected in ST theories since s = 1/2 for all BHs re- 
gardless of their masses or spins. For initially stationary 
scalar fields. Will and Zaglauer [7] proved that BBHs in 
Brans-Dicke theory are indistinguishable from binaries 
in GR to first post-Newtonian order. Similarly, Yunes, 
Pani and Cardoso [IS] recently extended this proof to 
all ST theories and to all post-Newtonian orders, but to 
leading order in the mass ratio. If the scalar field is not 
initially stationary, however, Horbatsch and Burgess [17] 
suggested that BBHs could retain scalar hair [TH] and 
emit dipole radiation, provided the holes have unequal 
masses. 

The goal of the present study is to investigate whether 
the conclusions from post-Newtonian and extreme-mass 
ratio studies regarding the BBH indistinguishability be- 
tween ST and GR theories carries over to the non-linear, 
comparable-mass regime. Our results confirm this indis- 
tinguishability, unless there exists a mechanism to induce 
dynamics in the scalar field. When the latter activates, 
the scalar field dynamics triggers dipole energy loss that 
leads to detectable differences in the and 'x' GW po- 
larizations. In the present study, we induce scalar field 
dynamics by placing the BBHs inside a scalar field bub- 
ble, which in some cases includes a potential. These in- 
homogeneities in the scalar field have a dramatic effect 
on the BBH dynamics and thus noticeable differences on 
the GW emission. We show that the changes on BBH dy- 
namics are due to accretion of scalar field by the merging 
BHs. 

We restrict attention to ST theories in vacuum with a 
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(Jordan frame) action of the form: 



S = 



2k 



(1) 

with K — 8tt G. Under a conformal transformation g^^^ — 
F{v)9piv^ or equivalently g^^ = {(f)) g , the action (in 
the Einstein frame) reads [31 13] : 

J d^x y/^ [R/{2 n) - V^</)V^</>/2 - y((/))] (2) 

where {d(j)/dip)^ = [(C/F) + (3/2)(F^/F)2] /k and V ^ 
U/{fiF^). We set ^(0) = e'^'^^'"^'/^ [31|3j. Thus, F = (/? 
and C — uj/(f with a; = —3/2 + K/{a — b(f>)'^. Brans-Dicke 
theory is recovered when 6 = 0, and GR when a = b = 0. 
We focus here on the case a — 0, and consider different 
values of b. 

The Einstein frame is convenient because the action 
in ([2]) yields the same equations as those of GR, namely 
G^^ = kT^u withT^^ = V^0V^0-(7^^ {Vp4>VP(l)/2 + V) 
and 00 = V^tj). Therefore, in the absence of a poten- 
tial, (j) — — constant yields G^i, — 0. Thus, vacuum 
spacetimes in ST theories will be equivalent to their cor- 
responding spacetimes in GR, independent of the choice 
of conformal factor ^4(0). This is also the case in the 
presence of a potential if one arranges for V{(j)Q) = and 
y,<t>{M = 0. We carried out simulations of spacetimes 
containing BH singularities that verify this exact equiv- 
alence between ST and GR, differing only at the level of 
round-off errors. 

Therefore, the only avenue to trigger differences be- 
tween ST theories and GR is by inducing dynamics 
in One possibility is with dtct) 7^ 0, e.g. 9t0 = 
et [17| . Another possibility is by introducing a direct 
scalar-field curvature coupling to the action, such as 
cj^R^^SpR^""^" mUS] or </>e"'^5pi?^.„/3i?'^''^'' [Il|20l|2l], 
that anchors the scalar field to the spacetime curvature. 
Yet another alternative is by introducing inhomogeneities 
in (/), which is the focus of the present study. 

We consider BBH inside a scalar field bubble with 
an initial profile 4'{r) — (/)o tanh [(r — ro)/cr] and an 
inflationary-inspired potential V = A(0^ — 0o)^/8 that 
yield discrete symmetry breaking |22| . The bubble has 
radius tq, thickness cr, and in the exterior (interior) 
4> = 4'o{~4'o)- Notice that, since we are interested in 
asymptotically flat spacetimes, (/)(r — 00) = (j)o', we shift 
the conformal factor as A = e^''^'^ ^'I'o)/^^ so the Einstein 
and Jordan frames are both Minkowskian and in the same 
coordinates ^3, ,4^ . This allows us to make fair compar- 
isons between BBHs in ST and GR theories, i.e. compar- 
isons among BBHs with the same BH masses, spins, sep- 
aration and eccentricity. If 7^ 0, the bubble becomes 
a topologically stable, domain wall with a « 2/{(j)oy/X) 
and surface energy density ^ « 3 0o-\/A/4 [22]. Since in 
the interior V{—(j)o) = and V,f,{—(l)o) = 0, the space- 
time in the neighborhood of the BBH can be arranged to 
be locally, nearly equivalent to its GR counterpart. 



Case 


00 


4 7rA/M2 Mo/M Mn/M 


a/Mh 


A 








0.990 0.952 


0.686 


B 


1/80 





1.179 0.963 


0.688 


C 


1/40 





1.747 0.999 


0.706 


D 


1/80 


10^ 


1.217 0.983 


0.685 



TABLE I: Parameters, masses and spin. 



We discuss results from a set of representative simula- 
tions in our study. The simulations were obtained with 
the Maya numerical relativity code of our group [231 121] j 
modified to include the scalar field in the Einstein 
frame. In all cases, the binary has non-spinning, equal- 
mass BHs in quasi-circular orbit, initially separated by 
11 M, with M the mass of the binary. The bubble sur- 
rounding the BBH has a radius rg — 120 M and thickness 
a = 8M. The simulations in the Einstein frame differ 
only in the parameters (/)o and A, as given in Table [T] 
Also in Table U are the values in the Einstein frame of 
the ADM mass Mq, the mass Mh and spin a/Mh of the 
final BH. Case A is the reference GR simulation, and D 
is the only case with non-vanishing potential. Extrac- 
tion of GWs is carried out in the Jordan frame. For each 
case, we extract GWs with values 6 = 0, 5 and 10 for the 
parameter in the conformal factor A = er^^'^ -<>o)/2_ 
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FIG. 1: Mass of the individual BHs as they accrete (j> before 
they merge. Here t = M is the initial time of the simulation. 

In all cases, the bubble shell collapses. For case D (a 
domain wall bubble), the thickness of the shell does not 
change much during its collapse. For cases B and G with 
I^ = 0, however, the shell disperses while collapsing. By 
the time the shell reaches the BBH, the cloud has en- 
compassed the binary. Independent of the details of the 
initial configuration, the scalar field has a dramatic effect 
on the BBH dynamics. When the bubble reaches the bi- 
nary, each BH accretes a significant amount of mass from 
the scalar field. Figure [l] shows the evolution of the mass 
(in the Einstein frame) of the individual BHs before they 
merge. The differences in the initial BH masses among 
all cases reflect the difficulty of setting up identical bi- 
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FIG. 2: Trajectories of the BHs when viewed from the Ein- 
stein frame. 



plotted in panels B, C and D correspond to 6 = 0. No- 
tice the obvious differences between the strains, and in 
particular when matched against the GR case A. How- 
ever, within each of the B, C and D cases, the differ- 
ences between a strain with b — and one with b ^ 
are undetectable. That is, the conformal factor does not 
add new features. We calculated mis-matches between 
the 5 = and 6 7^ strains for Adv. LIGO, within 
each B, C and D case, and found them to be < 10~'^. 
This can be understood by studying how the Weyl scalar 
^4 = —Ciffiifh conformally transforms. The tetrad vec- 
tors transform as v"- = A^^ v"" , and the Weyl tensor as 
Cabcd = Cabcd. and thus $4 = A-^ *4 = e^('^'-'^o) 
(tilde quantities correspond to the Jordan frame). Dur- 
ing the course of the simulation, we observed that the 
scalar field at the location of the GW extraction be- 
haves as (/) = 00 + '5</' with (5(/)/0o ^ 10^^. Therefore, 
^4 ~ (1 + 2 6(5<^0o)*4 « (1 + 10-^6 (/)g)*4- Since for 
the chosen parameters bcjx^^ 10^^, ^4 « ^"4 to one part 
in 10*. 



naries. Notice that for cases B and C the mass increases 
in two stages. The first increase occurs when the bubble 
is imploding and passes through the binary. The second 
increase is after the bubble bounces back and expands. 
The mass increase in C is larger because the bubble is 
more massive. With the choice of parameters, the most 
massive shell is C, followed by D and B (see the ADM 
mass Mo in Table]!]). For case D, the increase in the BH 
mass as a consequence of ^-accretion differs from cases B 
and C in the presence of multiple step-like increases. This 
is because in case D the bubble does not just bounce back 
and dissipate, but instead it lingers in the neighborhood 
of the binary, bouncing multiple times due to the poten- 
tial and leading to the series of step increases. Notice that 
for all cases, the BHs eventually reach a constant mass. 
The BHs gain 1%, 6% and 3% of their original mass for 
cases B, C and D, respectively. This gain in mass is cor- 
related with the mass Mh of the final BH. That is, in 
the GR case A, the final BH is Mh = 0.952 M, while in 
all other cases and to a good approximation, the final 
BH mass is Mh ~ 0.952 (M + 6M) with 5M the mass 
accreted by the merging BHs, i.e. 0.01 M for B, 0.06 Af 
for C, and 0.03 M for D (see e.g. Table]l|. 

The 0-accretion by the BHs has a dramatic effect on 
the binary dynamics. Figure ]2] shows the trajectories (in 
the Einstein frame) for all cases. Recall that case A is the 
GR reference case for a BBH in quasi-circular orbit. In 
all the other cases, the (/)-bubble induces eccentricity and 
accelerates the merger. The more massive the shell the 
larger the induced eccentricity. In particular, for C the 
influence of (j) is such that the binary basically plunges. 

Not surprisingly, the effect of the bubble is also felt 
in the emission of gravitational radiation. In Figure ]3] 
we plot the (2,2) mode of the GW strain polarization h+ 
(dashed line) and the amplitude \h\ — — ihy^ \ (solid 
line), with t = QM denoting merger time. The strains 
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FIG. 3; Mode (2,2) of the gravitational strain {r/M)h+ 
(dashed hne) and the amplitude {r/M)\h\ — {r /M)\h+ —ihx\ 
(solid Une). The time i = OM denotes merger time. 

In Figure ]3j cases B and D show the characteris- 
tic strain modulation observed in eccentric BBHs. The 
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waveform in case C, on the other hand, is basically a 
burst, since here the binary essentially plunges. Notice 
also that the amplitudes in C and D show small bumps, 
the first one at t ^ — 75 M for case C, and at t ^ —375 M 
for case D. These bumps are due to the shell as it goes 
through the binary. The bumps are therefore correlated 
with the mass jumps in Fig. [T] 
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FIG. 4: Mode (0,0) of the breathing mode, rM $22- In each 
panel, fe = 0, 5 and 10 are solid, dashed, and dotted lines, 
respectively, t — M denotes merger time. 

The presence of a scalar field in ST theories predicts 
dipole energy losses as well as a new GW polarization, 
sometimes called a breathing mode, given by the traceless 
Ricci tensor scalar, $22 — —Rimim = ~Rii/'^ [SH l26] . 
Figures |4] and [5] show respectively the modes (0,0) and 
(2,2) of $22 in the Jordan frame. In each panel, 6 = 0, 5 
and 10 are solid, dashed, and dotted lines, respectively. 
The time range in each panel was chosen to highlight 
when (f) interacts most strongly with the BHs, with t — 
M denoting merger time. The evident 6-dependence of 
$22 can be explained by how $22 conformally transforms: 
$2,2 = [$2,2 + bD{dt4i, ■■■)], with D a function that 
vanishes when dt(j) — 0. Therefore, at times when cf) un- 
dergoes significant evolution, $22 depends linearly on 6, 
as observed in Figs. ^ and [s] Overlaps of $22 with dif- 
ferent b in cases B and C correspond to moments when 
dt4> « at the wave extraction location. Then, as with 
\l/4, one has $22 ~ *i'22- For case D, this does not happen. 



i.e. <I>22 with different b values do not overlap because the 
potential induces longer lived dynamics in 0, persisting 
long after the binary has merged. 
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FIG. 5: Same as in Fig.|4]but for the mode (2,2). 



This study is a first step towards investigating de- 
tectable observational signatures in the GW emission 
from the late inspiral and merger of a BBH in ST gravity. 
Our results supports the view that, in order to "defeat" 
the no-hair constraint of BHs, and thus trigger detectable 
effects on the gravitational radiation, one needs a mecha- 
nism to excite scalar dynamics. In our study, we achieved 
this with inhomogeneities in the form of a scalar field 
bubble surrounding the BHs. The effect on the binary 
was dramatic for the chosen parameters, mostly due to 
the accretion of scalar field by the merging BHs. In a 
subsequent study, we will explore a broader range of pa- 
rameters, particularly those compatible with the param- 
eterized post-Newtonian bounds. We will also consider a 
situation with a cosmological scalar background, and the 
inspiral and merger of mixed binary systems. 
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